arXiv:1508.02208vl [physics.optics] 10 Aug 2015 


Resonant transmission of light in arrays of high-index dielectric nanoparticles 
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We study numerically, analytically and experimentally the resonant transmission of light in a 
waveguide formed by a periodic array of high-index dielectric nanoparticles with a side-coupled 
resonator. We demonstrate that a resonator with high enough Q-factor provides the conditions 
for the Fano-type interference allowing to control the resonant transmission of light. We suggest a 
practical realization of this resonant effect based on the quadrupole resonance of a dielectric particle 
and demonstrate it experimentally for ceramic spheres at microwaves. 
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Introduction. In the past years, a rapid development 
of the field of all-dielectric nanophotonics is observed 
aiming at the manipulation of optically-induced electric 
dipole (ED) and magnetic dipole (MD) resonances in di¬ 
electric nanostructures composed of nanoparticles with 
high refractive index m • Both ED and MD resonant 
responses were theoretically predicted mm and experi¬ 
mentally observed for silicon particles with linear 

size less than the operating wavelength, and it was shown 
that the resonance frequencies depend on the particles 
size and shape. Owing to these properties and to low 
losses in dielectric materials, different optical nanostruc¬ 
tures based on such nanoparticles were suggested M. 
One of such structures is a chain of dielectric nanopar¬ 
ticles with high refractive index. Such structures were 
suggested as a realization of optical waveguides with sub¬ 
wavelength light localization and low overall losses re¬ 
quired for a design of highly efficient integrated photonic 
circuits 0E1- 

A common practical scheme that employs a coupling 
between an optical waveguide and a high-Q cavity was 
suggested for controlling the light propagation and real¬ 
ization an optical bistable devices jT2FFT8] . Most of the 
studies were carried out for photonic-crystal waveguides 
with different geometries and different types of defects or 
cavities m mi qshu]- In such a case the interference 
between the continuum of states (a transmission band of 
the waveguide) and a discrete resonant state (a resonance 
of the cavity) provides the conditions for the Fano-type 
resonance [23] realized in a sharp and asymmetric line 
shape, which reduces the frequency shift and the value 
of energies required for all-optical nonlinear switching, 
logic, and information processing [24, 25.• 

In this Letter, we study the resonant transmission of 
light through a waveguide formed by a chain of high- 
index dielectric nanoparticles with a side-coupled defect 
nanoparticle. Localization in such structures in trans¬ 
verse direction is due to the high index of dielectric ma¬ 
terial and its guiding properties originate from the long- 
range coupling between the particles. Subwavelength size 
of the waveguide and the possibility of controlling light 


transmission makes it promising for designing the effi¬ 
cient optical integrated circuits. The results of this work 
can be applied to any spectral range, if appropriate ma¬ 
terials are provided. Following this idea, we conduct mi¬ 
crowave experiments with ceramic spheres and confirm 
all predicted resonant phenomena. 

Theoretical model. We consider a chain of closely 
spaced high-index dielectric nanoparticles with side- 
coupled defect nanoparticle as shown in Fig.[l][a). We fix 
the permittivity of dielectric particles Re(e p ) = 16 with 
small material losses Im(e p ) = 0.02, that approximately 
corresponds to silicon nanoparticles in optics |26| and to 
ceramic spheres at microwaves m , and the period of the 
chain a = 2 R. The radius Ry , and the permittivity Ed of 



Figure 1: (a) Schematic of a chain of high-index dielectric 
nanoparticles with permittivity e v , radius R and period a 
coupled to a side particle with permittivity £d and radius Rd. 
“A” is a position of the source and “B” is the position of the 
detector, (b) Dispersion curve for the transversely polarized 
mode of an infinite chain of spherical nanoparticles; /3 is the 
Bloch wavenumber, A is the wavelength. 
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the side-coupled particle vary for different cases. 

For theoretical description of the considered system we 
use analytical approach when applicable, and direct nu¬ 
merical calculations in other cases. Analytical approach 
is based on the well-known coupled-dipole model, where 
each particle is modeled as a hybrid of magnetic and elec¬ 
tric dipoles with magnetic m and electric p momenta, 
oscillating with frequency uj [oc exp(-iut)]. In the CGS 
system of units such approach formulates as follows: 

Pi = Olei [pi]Pi ~ , 

m, = a mi Y (d) m i + . 

3 

where a m and a e are the magnetic and electric polar¬ 
izabilities of a sphere [2Z], Cij — AijI + (8) Yij), 

Gij = —DijYij x /, (8) is the dyadic product, I is the unit 
3x3 tensor, r^- is the unit vector in the direction from 
i-th to j -th sphere, and 

A _ exp (ikhVij) f l2 1 , ik h \ 

u _ exp (ik h rij) ( ;2 , 3 3 ik h \ 

ij ~ r tj { h + (rij ) 2 rij ) ’ 

D = exp (ikhrjj) / fc2 iJY \ 

13 n, V h r^ J ’ 

where is the distance between the centers of i-th and 
j-th spheres, is the permittivity of the host medium, 
kh = v / ^ cj / c the host wavenumber, uj = ‘htv, v is 
the frequency, and c is the speed of light. It is known, 
that such model provides very accurate results for wide 
range of parameters m In the cases when the dipole 
model becomes inapplicable, numerical calculations are 
performed in CST Microwave Studio. 

Numerical results. We operate in the frequency range 
near the magnetic dipole resonance frequency (which is 
about 2R/X « 1/y / Re(e p ) = 0.25) where magnetic mo¬ 
ments are dominant. In Fig. [ljb) we show the disper¬ 
sion curves for the transversely (i.e. when magnetic 
dipoles are oscillating in direction perpendicular to the 
chain axis) and longitudinally polarized (when magnetic 
dipoles are oscillating in direction parallel to the chain 
axis) magnetic modes of the infinite chain with the cho¬ 
sen parameters (assuming e p to be real-valued). Further 
calculations are carried out for the TM modes. 

First, we use the side-coupled particle with permittiv¬ 
ity £d = £p = 16 + 0.02i. The MD resonance frequency of 
the particles with the chosen parameters lies outside the 
TM band of the uniform chain [see Fig. [ljb)] . By increas¬ 
ing the radius of the side particle to Rd = 1.1 i?, we shift 
its resonance frequency uJd to the transmission band. We 
excite the chain with a point magnetic dipole placed at 
the distance of 0.5 R from the left end of the chain [point 
“A” in Fig. [lja)] and calculate the intensity of the i-th 
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Figure 2: Transmission through the chain of 60 particles for 
different parameters of the side-coupled sphere: Ed — £ P , Rd — 
1.1 R [dashed red curve (2) - z-excitation; dashed blue curve 
(3), y-excitation]; Ed — 64, Rd — 0.551? [solid red curve (4) - 
^-excitation; solid blue curve (5), y-excitation]. Black dashed 
curve (1) corresponds to the case of the chain without side- 
coupled particle. 


component of the magnetic field \Hi\ 2 at the other side of 
the chain [at the distance 0.5i? from the right end, point 
“B” in Fig. [lja)]. Which field component is detected is 
defined by the polarization of the excited chain modes. 
In this geometry two polarizations can be distinguished: 
(m x ,m y ,p z ) and (px^Py^z)- The first type of modes 
can be excited with the ^-oriented magnetic dipole and 
therefore we refer to this case as to the z-excitation. We 
excite the second type of modes with ^/-oriented mag¬ 
netic dipole (^/-excitation), since longitudinal x-polarized 
modes lie in different spectral region and we do not con¬ 
sider them here. In this view the case of ^-excitation is 
the simplest, since magnetic moments in all particles are 
oriented in ^-direction. But in the case of the ^/-excitation 
magnetic dipoles oscillate in x — y plane and the coupling 
between the particles is stronger, therefore the stronger 
effect may be observed. 

In Fig. [ 2 ] we show \Hi\ 2 at the point “B” for the chain 
with and without side particle. Here “i” is y or z for y— 
or z-excitation, respectively. The side particle is located 
at the distance yd = 2.5 R from the chain. In the case 
of the ^-excitation (dashed red curve (2) in Fig. [ 2 ]) the 
transmission remains almost the same, and in the case of 
the //-excitation (dashed blue curve (3) in Fig. [ 2 ]), when 
the coupling is stronger, we observe a very broad dip that 
is stretched over the whole pass band. From these calcu¬ 
lations it is clear that the Q-factor of a single particle is 
too low to provide some sharp features that can be ob¬ 
served, neither for the z-excitation, nor for //-excitation. 
In order to show that side-coupled particle can in fact 
provide the conditions for Fano-type resonance, we in- 
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Figure 3: Relative transmission for through a chain with side- 
located sphere with Ed — 64 and Rd — 0.55R; curves (1) 
correspond to the case ^-excitation, curves (2) to the case of 
//-excitation. Dashed red curves are obtained via the calcula¬ 
tions made within a dipole model and solid black curves via 
formula 0. 


crease its Q-factor in the simplest theoretical way - by 
increasing its permittivity to Sd = 64. We also decrease 
the radius of side-coupled particle to Rd = 0.55i?, so 
its MD resonance frequency remains approximately the 
same. In this case, as it is shown with solid red curve 
(4) and blue curve (5) (for z- and //-excitation, respec¬ 
tively) in Fig. [2} we observe the sharp resonance dip in 
the transmission spectrum of the chain. 

Further we investigate the case of high-Q side-coupled 
particle with Ed = 64 in more details. Transmission 
coefficient through a single-mode waveguide with sym¬ 
metrically side-coupled resonator, that supports n eigen- 
modes with the resonance frequencies Q n and the mode 
decay rate can be calculated with the following for¬ 
mula P3M3S]: 


* = i-£ 




Au n + i(r° + n n ) 


; T = \t\\ (l) 


where Acj n = uj — fi n , and is the decay rate of the 
n-th mode of the side-coupled resonator into the waveg¬ 
uide. Here relative transmission is determined as the ra¬ 
tio of the magnetic field intensities at the detector point 
T = \Hf \ 2 /\Hi\ 2 with and without side particle. Note 
that the side particle also creates the field in the detector 
point directly through the free space, however this contri¬ 
bution is noticeable only for short chains < 20 particles. 
Our calculations are made for the chains consisted of 60 
particles, therefore this contribution is negligibly small. 

In Fig. [3] show the relative transmission T for the cases 
of z- and y- excitation. Solid black curves are calcu¬ 
lated with dipole model, and dashed red curves with for¬ 
mula 0. In the case of z- excitation waveguide mode 
induces only ^-oriented magnetic moment on the side- 


coupled particle, i.e. it interacts only with one MD eigen¬ 
mode, and transmission coefficient reduces to m 


Au 2 + (r°) 

Aw 2 + (r c + r 0 ) 2 ’ 


In this case almost symmetric dip in transmission is very 
well described by the formula § for the picked parame¬ 
ter r c . 

In the case of y- excitation the waveguide mode can 
propagate via excitation of both x- and //-oriented mag¬ 
netic moments on the side-coupled particle. In this case 
the coupling is stronger and the dip in transmission is 
deeper. We also observe that the dashed symmetric curve 
(2) in Fig. [ 3 ] fitted with formula ([l]) can not describe the 
transmission very accurately, since the solid curve (2) 
is assymetric. For more accurate desription the depen¬ 
dence of T c on the frequency, that causes the assymetry, 
has to be taken into account. Calculations presented in 
Fig. I confirm that the resonance features in transmis¬ 
sion spectra appear as a result of Fano-type interference 
of the waves that go through resonant and non-resonant 
channels of transmission. 

Experimental results . In order to verify this concept ex¬ 
perimentally, we scale the dimensions to the microwave 
frequency range and perform an experimental study of 
the transmission efficiency of the waveguide composed 
of 30 identical ceramic disks with radius R = 3.95 mm 
and height h = 3.9 mm and permittivity e ~ 16.4. The 
disks are placed in a holder made of a styrofoam mate¬ 
rial with permittivity close to 1. The axis of the chain 
is oriented perpendicular to the axis of the disks. The 
magnetic dipoles are realized as copper loops with the 
diameter of 2 mm being connected to Agilent E8362C 
Vector Network Analyzer via coaxial cables. Transmit¬ 
ting and receiving loops are placed at the distance of 1 
mm above the first and last disks, respectively. 

A practical realization of the resonator with a high 
value of the Q -factor at the subwavelength scale is a 
complicated task mostly due to radiation losses. As a 
realistic case, here we suggest to use a sphere with larger 
radius that exhibits higher-order multipole resonances, 
with the frequencies tuned to the operational frequency 
by changing its size. It is known that a dielectric sphere 
with e ~ 16 and radius R exhibits a magnetic quadrupole 
(MQ) resonance at the frequency ~ 0.175 c/R, with the 
Q-factor being substantially higher than that of MD and 
ED resonances [4]. For a sphere with Rd = 6.1 mm MQ 
resonance is expected at the frequency « 8.6 GHz. In 
our experiments the permittivity of the spheres was a bit 
smaller than 16, therefore the MQ resonance frequency 
is slightly higher. 

The magnitude of the experimentally measured trans¬ 
mission coefficient (|S f 2 i| parameter) is shown in Fig.^a). 
Black and red curves correspond to a chain without and 
with side-coupled sphere, respectively. The dip in the 
transmission (a red curve) is observed at the frequency 
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Figure 4: (a) Experimentally measured and (b) numerically calculated transmission efficiency for the chains of ceramic disks 
without (black curve) and with (red curve) side-coupled sphere, (c-e) Magnetic field distribution \H Z \ in the chain with 
side-coupled sphere with Rd — 6.1mm for three frequencies: (b) v — 8.5 GHz, (c) v — 8.725 GHz, (d) v — 9.0 GHz. 


of the MQ resonance of the side-located sphere. For other 
frequencies, the transmission efficiency remains approxi¬ 
mately the same. We support these experimental results 
by direct numerical simulations performed with the help 
of the CST Microwave Studio, since in this case the an¬ 
alytical dipole model is no longer applicable. The cal¬ 
culated transmission for the same parameters as in ex¬ 
periment is presented in Fig.[4^b). We observe very good 
qualitative agreement between numerical and experimen¬ 
tal results. A quantitative difference between the exper¬ 
imentally measured and numerically calculated spectra 
can be explained by a noticeable dispersion of the size 
and permittivity of real ceramic disks and by difficulties 
of an accurate modeling of the transmitting and receiving 
loops. 

In Figs. [4^ c-e) we show the calculated magnetic field 
distribution \H Z \ in the plane y = 0 at three different fre¬ 
quencies. We observe that an additional particle almost 
completely blocks the propagating wave, which in its turn 
induces a multipole moment in the particle. Since the 
considered system is open, a weak scattering of radiation 
is also observed. 

In summary, we have demonstrated theoretically, nu¬ 
merically and experimentally the possibility of an res¬ 
onant control of the light transmission through waveg¬ 
uides composed of high-index dielectric nanoparticles 
with side-coupled resonators. We have demonstrated 
that an effective control requires a resonator with high 
enough Q-factor, that can be realized for a dielectric 
sphere with a larger radius supporting high-Q multipole 
resonance in the operating spectral range. 
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